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Abstract 

The electrical resistance of ferromagnetic/normal-metal (F/N) heterostructures depends on the 
nature of the junctions which may be tunnel barriers, point contacts, or intermetallic interfaces. 
For all junction types, the resistance of disordered F/N/F perpendicular spin valves as a function 
of the angle between magnetization vectors is shown to obey a simple universal law. The spin- 
current induced magnetization torque can be measured by the angular magnetoresistance of these 
spin valves. The results are generalized to arbitrary magnetoelectronic circuits. 



PACS numbers: 72.25.Ba,75.70.Pa,75.60.Jk,75.75.+a 



Magnetoelectronics achieves new functionalities by incorporating ferromagnetic materials 
into electronic circuits. The giant magnetoresistance, i.e. the dependence of the electrical 
resistance on the relative orientation of the magnetizations of two ferromagnets in a ferro- 
magnetic/normal/ferromagnetic (F/N/F) metal structure or "spin valve", is applied in read 
heads of high information density magnetic storage systems Usually, such a device is 
viewed as a single bit, the magnetizations vectors being either parallel or antiparallel. Early 
seminal contributions by Slonczewski || and Berger || revealed fundamentally new physics 
and technological possibilities of noncollinearity, which triggered a large number of experi- 
mental and theoretical research. An important example is the non-equilibrium spin-current 
induced torque (briefly, spin torque) which one ferromagnet can exert on the magnetization 
vector of a second magnet through a normal metal. This torque can be large enough to 
dynamically turn magnetizations ||, which is potentially interesting as a low-power switch- 
ing mechanism for magnetic random access memories ||. The spin torque is also essential 
for novel magnetic devices like the spin-flip transistor || [7|, ||, detection of spin-precession 



, the Gilbert damping of the magnetization dynamics in thin magnetic films [10|, and 



spin-injection induced by ferromagnetic resonance [1 I 



Recently, two theoretical approaches have been developed which address charge and spin 
transport in diffusive noncollinear magnetic hybrid structures. The magnetoelectronic "cir- 
cuit theory" |5j is based on the division of the system into discrete resistive elements over 
which the applied potential drops, and low-resistance nodes at quasi-equilibrium (as in 
Fig. 1(a)). The electrical properties are then governed by generalized Kirchhoff rules in 
Pauli spin space and can be computed easily. Each resistor is thereby characterized by four 
material parameters, the spin-up and spin-down conductances = ^2 nm [S nm — |^lm | 2 ] as 
known from the scattering theory of transport |Ll|, as well as the real and imaginary part 
of the "mixing conductance" = ^2 nm [S n m — r l m ( r nm)*]> where r s nm is the reflection coef- 
ficient between n— th and m— th transverse modes of an electron with spin s in the normal 



metal at the contact to a ferromagnet. Waintal et al. studied the random matrix theory 
of transport in noncollinear magnetic systems as sketched in Fig. 1(b). Their formalism 
did not require the assumption of highly resistive elements, but the algebra of the 4x4 
scattering matrices in spin space seemed so complex, that analytical results were obtained 
in limiting cases only. 

Both theories are not valid in the limit of intermetallic interfaces in a diffuse enviro- 



2 



ment (see Fig. lc) like the perpendicular spin valves, studied thoroughly by the Michigan 
State University collaboration |TJ]] and others []I5 , [T6|j . These studies provided a large body 



of evidence for the two-channel (i.e. spin-up and spin-down) series resistor model and a 
wealth of accurate transport parameters like the interface resistances for various material 
combinations. Transport through transparent interfaces in a diffuse environment has been 
studied for collinear magnetizations by Schep et al. W7\ . Under the condition of isotropy of 



scattering by disorder, it was found that the bulk resistances, which are proportional to the 
layer thicknesses, are in series with interface resistances, for each spin s 

~g s ~ g s 2 \N* + N N J ' U 

where Nf and are the number of modes of the bulk materials on both sides of the F/N 
contact. Physically, in Eq. ([I]) a spurious Sharvin resistance is substracted from the result of 
scattering theory. This correction is large for transparent interfaces and essential to obtain 
agreement between experimental results and first-principles calculations |T7], [15, [19 . 



In exchange-biased spin valves, it is possible to measure the electric resistance as a func- 
tion of the angle between magnetizations, which has been analyzed experimentally and 
theoretically pD , |22"f . The present study has been motivated by Pratt's observation that 



this angular magnetoresistance could accurately be fitted by the form 

R(9)-R(0) _ 1-cosfl 
R(n) -R(0) ~ x(l + cos#) + 2 

with one free parameter x that is given by circuit theory 



(2) 



» = (3) 
1 — p z tier] 

in terms of the normalized mixing conductance 77 = 2g^/g, the polarization p = (g^ — g{) /g, 
and the average conductance g = g^ + g^. This was surprising, since the circuit theory, as 
mentioned above, was not designed for metallic multilayers, and, indeed, the numerical value 
of fitted parameters did not make sense, also after including effects of bulk scattering in the 
ferromagnetic layers ||23|| . 

In the following we develop a theory of transport in disordered magnetoelectronic circuits 
and devices in the diffuse regime which unifies and extends previous theoretical approaches. 
We find simple analytical results with parameters that are accessible to realistic electronic- 
structure calculations. The angular magnetoresistance for perpendicular spin valves agrees 



with the universal form [Eq. @] in agreement with measurements [22]|, and is used to 



determine the mixing conductance and spin torque. The theory is valid under two conditions: 
(i) the system should be diffusive, i.e. the elastic mean free path £ (including scattering at 
interfaces) should be smaller than typical sample scales and (ii) the ferromagnetic elements 
should have an exchange splitting A, which is large enough that the magnetic coherence 
length l c = h/y/2rnA < min (£, (If) , where <1f is the thickness of the ferromagnetic layer. 
These conditions are usually fulfilled in transition-metal systems: Deviations from diffusive 
behavior, like quantum-size effects and breakdown of the series resistor model, are small 



or controversial |24fl , whereas the magnetic coherence length is of the same order as 



the lattice constant in high-T c transition-metal ferromagnets [10, 25]. We obtain identical 



results by two methods: The first one is a combination of the Boltzmann-like method of 



Schep et al. [17] f° r collinear systems and the random-matrix theory of Waintal et al. [|13 
The second one is an extension of magnetoelectronic circuit theory |J to arbitrary resistors. 

Let us consider planar spin- valve structures as shown in Fig. 1. We assume the existence 
of a distribution function at a certain position x in the sample (a "node"), which in spin- 
polarized systems has 8 elements f ss , (x). We arrange them into a 4 x 1 vector f ± = 
(/rp ft\i ft\i ft\) T as we U as a 2 x 2 matrix, denoted by a hat: 

f± , x) = I (*) (*) ] (4) 

V At (*) fti (*) J 

The superscript denotes that the distribution is in general anisotropic in reciprocal space, + 
for right-moving — for left moving, indicating that, in contrast to || the current density 
in the nodes is not negligible. The distribution functions at different nodes are matched via 
boundary conditions: 

/+ (x B ) = f A ^ B f + (x A ) + Rb-*b?~ (x b ) (5a) 
f- (x A ) = R A ^ A f + (x A ) + Tb^aT M , (5b) 

where the 4x4 transmission and reflection probability matrices (indicated by the caret) 
have elements like [ |l^| : 

P^-Jy = J^bYI (^nm 5 )] (6) 

1 nm 

where = (S^i + 8^2) + (6^3 + 8^) , iVf is the number of modes for spin s in B, 
and tnm B * s a vector of the transmission coefficients in spin space. 



Let us calculate the electrical charge current in a symmetric two-terminal spin valve with 
relative magnetization angle 9 (Fig. 1). xl and xr are within left and right ferromagnets 
at a distance from the interface equal to the spin diffusion length in the ferromagnet £f d 3> 
£ c , and thus define the magnetically active region. In the coordinate systems defined by 
the magnetization directions, the transverse components of the spin accumulation in the 
ferromagnets vanish [0, and the distributions in the magnets depend on the local spin 
current densities 7 S and (spin-independent) chemical potentials \x only: 

/± (x) = ((± 7t + /i) (x) , 0, 0, (±7 A + n) (x)) . (7) 

In symmetric junctions the spin current is symmetric as well, 7 S (xl) = 7s (xr). The charge 
current i c = (e 2 /h) J2 S N^ 7 S divided by the chemical potential drop equals the electrical 
conductance G = i c /Afi. Eqs. (|5],|7|) then lead to: 

G = yEKp- f L^R + Rr^r] ^ T L -> R } i . • (8) 



=1,4 



In principle, the matrices T and R do not need to be approximate. 

In dirty systems, more nodes may be introduced at convenient locations in the sample and 
Eqs. (|5p implies that total transport probability matrices can be composed in terms of those 
of individual elements by semiclassical concatenation rules [BBI]. For instance, the transmis- 
sion through a F(0)/N/F(#) double heterojunction as in Fig. 1 (without bulk scattering) 
takes the form: 

f (6) = T N ^ F (6) [1 - R N ^ N (0) i? N ^ N (6)} - 1 T F ^ N (0) . (9) 

These rules have been derived from the (phase-coherent) scattering theory by averaging over 
random matrices [13|] and found to be valid to leading order in iV^ 1 , where A^v is the number 



of transport channels in the normal metal. Bulk impurity scattering can be represented by 



diagonal matrices |L3], [17 



1 e 2 pfd 



N? h A 



B 



-1 



where pf , ds, Ab are the bulk resistivities, thickness and cross section of the bulk material 
B, repectively. 
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The problem can be simplified by transformations into the coordinate systems defined by 
the magnetization directions of the ferromagnets. In terms of the spin-rotation 

(cos | — sin | \ 
sin § cos | J 

and projection matrices (s = ±1) 

1 / 1 + s cos 6 s sin 9 \ 

Us(e) = -[ , (12) 

z \ s sin 6 1 — s cos 6 I 

the interface scattering matrices (omitting the mode indices for simplicity) are transformed 
as i t%7 N = U ss df, tf s r F = tfUl, r N ^ N = £> s r s cJV , and r F s ~T F = rf5 ss >, where the 
superscript c indicates that the matrices should be evaluated in the reference frame of the 
local magnetization and spin-flip scattering in the contacts has been disregarded. 

The angular magnetoresistance can now be evaluated analytically for our spin valve in 
terms of the three interface conductances g^, g^, defined above, the bulk number of 
modes Nf, N , and bulk resistances p F , p , whereas the magnetization angle and layer 
thicknesses are the variables. Surprisingly, the form Eq. (|2|) is recovered, but with renor- 
malized parameters. The spin-dependent interfaces conductances are identical to Eq. ([[]), 
whereas, including the effect of bulk scattering, 

9n 9n 2 \ A n n) ' 

By letting Nf — > oo we are in the regime of ||13|| . The circuit theory is recovered when, 



additionally, Nn — > oo. The bare mixing conductance is bounded not only from below 
R> e 9ll ^ 9/2 0) but also from above \g-\i\ 2 /Rep^ ^ 2N N - The polarization and relative 
mixing conductances are also renormalized, with < \rj\ < 00. 

It is not obvious how these results should be generalized to more complicated circuits 
and devices and to the presence of spin-flip scattering in the normal metal. The mag- 
netoelectronic circuit theory does not suffer from these drawbacks. In the following 
we demonstrate that above results can be obtained with less effort, proving that with the 
renormalization of the transport parameters by subtracting Sharvin resistances, circuit the- 
ory remains valid for arbitrary contacts. To this end we construct the fictitious circuit 
depicted in Fig. 2. Consider a junction which in conventional circuit theory is characterized 



6 



by a matrix conductance g leading to a matrix current % when the normal and ferromagnetic 
distributions fi and fn are not equal. When the distributions of the nodes are isotropic, we 
know from circuit theory that: 

^^(^ss'Us^h- fn)u sl , (14) 



fs,—s- 



where the projection matrices it s are defined in Eq. ( |T2"D and (g) ss = g s , (g) s ,-s =9s 
Introducing lead conductances, which modify the distributions }l — > fi and f 2 <— fa, 
respectively, we may define a (renormalized) conductance matrix g, which causes an identical 
current i for the reduced (matrix) potential drop: 

i = ^2 (f) (A - /a) u s > • (15) 

ss' 

When the lead conductances are now chosen to be one-half of the Sharvin conductances, 
and using (matrix) current conservation: 



2N N \J L - A J (16) 
522N?u. (} 2 - f R ) (17) 



straightforward matrix algebra leads to the result that g is identical to the renormalized 
interface conductances found above [Eqs. (pO) and, without the bulk term, (O)]. By replacing 
g by g we not only recover results for the spin valve obtained above, but we can now use 
the renormalized parameters also for circuits with arbitrary complexity and transparency of 
the contacts. Also spin- flip scattering in N can be included |J; it does not affect the form 
of Eq. @) either, but only reduces the parameter x- 

Experimental values for the parameters for Cu/Permalloy (Py) spin valves are x = 1-2 
and p = 0.6 p2[ . Disregarding a very small imaginary component of the mixing conductance 
0, using the known values for the bulk resistivities, the theoretical Sharvin resistance for Cu 
(0.55 • 10 15 fi _1 m _2 /spin [fL7|l), and the spin- flip length of Py as the effective thickness of the 
ferromagnet Uf d = 5 nm we arrive at the bare Cu/Py interface mixing conductance 

G n = 0.39 (3) ■ 10 15 Vl^m" 2 , which is close to that of Co/Cu §. 

The analytical expression for the spin torque on either ferromagnet, i.e. the spin current 
normal to the magnetization direction, reads 

L(6) = P -l 7 ^ ^, (18) 

v ; 2 (77- l)cosfl + l + ?7 2tt y j 
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in terms of parameters which can be measured as well as computed from first principles. 



Previous results H, 13] are recovered in the limit that r\ — > 2 and p — > 1. By the generalized 



circuit theory it is straightforward to compute the torque on the base contact of the spin- 
flip transistor with antiparallel source-drain magnetizations (three identical contacts) ||. 
Interestingly, it is larger and has a symmetric and flatter dependence on the angle of the 
base magnetization direction 9 : 

= pgr) sin 6 A/x 

61 ' (l-fj)cos 2 e + 2 + fi 2n ' 

This opens the way to engineer materials and device configurations to optimize switching 
properties of magnetic random access memories. 
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FIG. 1: Different realizations of perpendicular spin valves, (a) Highly resistive junctions like point 
contacts and tunneling barriers limit the conductance, (b) Spin valve in a geometrical constric- 
tion amenable to the scattering theory of transport, (c) Magnetic multilayers with transparent 
interfaces. 9 is the angle between magnetization directions. 

FIG. 2: Fictitious device which illustrates the generalization of circuit theory to transparent resis- 
tors as discussed in the text. 
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